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Abstract 

We propose a mechanism displaying gluon confinement, as denned by the behavior of the propaga- 
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theory giving rise to a replica model composed of three mixed 577 (2) groups. The mechanism consists 
in the usual 577(3) Yang-Mills theory in the Landau gauge, with a soft breaking term in such a way as 
to change the field propagation and group content at low energies. The relation of this soft mass term 
with the Gribov problem is presented and the link between soft terms and the scaling and decoupling 
solutions is discussed. 
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1 Introduction 



The problem of gluon confinement is the core of the general investigations of strongly coupled gauge 
theories. An important aspect of gluon confinement is related to the behavior of the gluon propagator. 

Recently this problem has received great attention in different approaches. One of these approaches 
comes from lattice simulation where the behavior of the gluon propagator in the infrared regime is studied 
[T| 1281 13]. These results display positivity violation thus making impossible a particle interpretation for 
the gluon excitation at low energies. This is taken as a strong signal of gluon confinement. In the 
analytical point of view, one possible approach of the confinement problem comes from the analysis of the 
Gribov copies |4j, where the Gribov-Zwanziger (GZ) model and this refined version, the so-called 

Refined Gribov-Zwanziger (RGZ) model [7J, take place. Also, a recently developed model based on the 
introduction of a replica of the Faddeev-Popov action enjoys a confined gluon propag atoiH 181. Usually, 
these models provide propagators behaving a^l 

where D(p 2 ) is the gluon form factor in Euclidean spacetime and 7 is a mass parameteid. It is easy to 
notice that a propagator of this type has complex poles, being impossible its identification with a propa- 
gation of a physical particle. In other words, it is a suitable candidate to be a confined object. 

An important feature that has been studied about propagators like (UJ) is that they can be seen as a 
propagation of two unphysical modes with imaginary squared masses iry 2 , named z-particles [9]. In fact, 
one can immediately notice that the form factor ([I]) can be written as 

^ =u^— 2+ * y (2 ) 

p A + 7 4 2 \p 2 - *7 P 2 + H J 

Despite the fact that such propagators do not have an interpretation in the physical spectrum, it is still 
possible to construct composite operators, whose the correlation functions exhibit a Kallen-Lehmann 

spectral representation: 

(0(p)0(-p))= [°°dT-?^L, (3) 



TO 



T + p 



where p(r) is the positive spectral density and tq > stands for the threshold. An important feature of 
expression ([3]) is that we can move from Euclidean to Minkowski space. Moreover, the positivity of the 
spectral density p(r) enables us to give an interpretation of ([3]) in terms of physical states with positive 
norm. 



We will refer to this model just as the replica model. 

2 In the case of the RGZ model the gluon propagator take into account the effects of dimension two condensates and it is 
more complex. Typically, the gluon form factor looks like 

D( 2, = p 2 + M 2 

(P ' p 4 + (M 2 + m 2 )p 2 + ( 7 4 + A'Pm 2 ) ' 

where M and m are mass parameters associated with the condensates, while 7 is the so-called Gribov mass parameter. 

3 In the GZ model this parameter is known as the Gribov parameter, which is directly associated with the restriction of 
the Feynman path integrals to the Gribov region. 
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In the present work, we will study a SU(3) Euclidean Yang-Mills theory and try to explore a well known 
property of this group, which consists in the fact that it contains three mixed SU(2) groups [TOj [11] . In 
the presence of a dimension two condensate, {d abc A^A^) , with d abc being the totally symmetric structure 
constant of the SU(3) Lie algebra, the color group symmetry is explicitly broken. However, two of these 
three ST/ (2) groups are related in a structure of i-particles, and, thanks to a residual symmetry, it is 
possible to write composite operators having a spectral representation like ([3]). 

The paper is organized as follows. In Sect. 2, we make a brief review about the GZ model, the in- 
troduction of the z-particles, and the replica model. In Sect. 3, we explain how the three mixed SU(2) 
groups can emerge from the SU(3). In Sect. 4, we introduce the condensate (d abc A b A c ^ and write the 
action corresponding to the model that we desire to study. Also in this section we show how this con- 
densate modifies the gluon propagator and how the z-particle structure appears between two of the three 
SU(2) groups. In Sect. 5, the relation of the mass terms, obtained from the introduction of dimension two 
condensates, and the scaling and decoupling solution is established. The Sect. 6 is dedicated to discuss the 
Kallen-Lehmann spectral representation of the candidates to be a physical observable. Our conclusions 
are presented in Sect. 7. The symmetry content of the model, characterized by a full set of Ward identities 
compatible with the Quantum Principle Action [12] , and the proof of its renormalizability are presented 
in the Appendix A. 

2 A brief review: the GZ model, i-particles and the replica model 
2.1 The GZ model 

In order to clarify the understanding of the Gribov problem and its correlation with our proposal we will 
present below a description of the problem in the Landau gauge. The Euclidean SU(N) Yang-Mills action 
in the Landau gauge is given by: 

S YM = J d 4 x (jF^F- v + ib a d^A- + cfd,D ab c h ^ , (4) 

where 

F; v = d,A a u - d„A« + gr bc A\Al (5) 

and 

Df = S a %-gf abc A*. (6) 

Here, A^ is the gauge field, b a is a Lagrange multiplier enforcing the Landau gauge, d^A^ = 0, (c a ,c a ) 
are a pair of anti-commuting scalar fields known as the Faddeev-Popov ghost fields, and g is the coupling 
constant of the theory. The labels (a, b, c, . . .) run to 1 to (A 2 — 1) and f abc are the totally anti-symmetric 
structure constant of the Lie algebra of the generator of SU(N). Also, This action is left invariant under 
the following nilpotent BRST transformations: 

sAl = -D ab c b , sc a = ^f abc c b c c , s^ = ib a , sb a = 0. (7) 

Although the gauge be fixed by the Faddeev-Popov method, Gribov showed in [4"pl that there are still field 
configurations obeying the Landau gauge linked by gauge transformations, i.e. there are still equivalent 
configurations, or copies, being taken into account into the Feynman path integral. In other words, the 

4 See also [13] for a pedagogical review. 
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gauge is not completely fixed and the remaining ambiguity is allowed due to the existence of normalizable 
zero-modes of the Faddeev-Popov operator, 

M ab = -d^Df . (8) 

Gribov also showed that to eliminate these copies the domain of integration of the functional integral 
should be restricted to a certain region f2, the so-called Gribov region, that is defined as the set of field 
configurations performing the Landau gauge condition, for which the Faddeev-Popov operator is strictly 
positive, namely 

n-.= {A-\d fl A^ = 0,M ab (A)>0}. (9) 

Its boundary, d£l, where the first vanishing eigenvalue of the Faddeev-Popov operator shows up, is known 
as the Gribov horizon. 

As in the region f2 the Faddeev-Popov operator is positive than its inverse must diverge when approach- 
ing the horizon, due to the existence of a zero mode. So the restriction to the first Gribov region is 
implemented requiring that 

' " " ->afc\-l| 



G(p, A) = , (10) 

which is the normalized trace of the ghost connected two point function in momentum space, has no pole 
for a given nonvanishing value of the momentum p, except for the singularity at p = 0, corresponding to 
the first Gribov horizon. At p ~ one can write 

1 1 



G(p,A) 



p 2 1 — a(p, A) 



From the above expression (lllj) . it follows that the no-pole condition at finite nonvanishing p is 

a(p,A)<l. (12) 
As a(p, A) decreases as p 2 increases one can also take 

It is important to emphasize here that we work with the trace of the ghost propagator to find the restric- 
tion to the Gribov region. This is a particularity of the Gribov mechanism in the Landau gauge. This 
is related to the convexity of the Gribov region in the Landau gauge. Other gauges, like the maximal 
Abelian gauge, does not necessarily present the same property [HI CE21 [160. 

In order to perform the restriction to the Gribov region into the partition function, Z, the final step 
is to introduce the no-pole condition with the help of a Heaviside function: 

Z = j VA5(dA)6(l - <r(0, A)) exp~ SYM . (14) 

This will give rise to a propagator for the gauge field of the type 

W-P)4(P)> = (V " V ~f) , (15) 



3 In the maximal Abelian gauge we take only the trace of diagonal ghost propagator. 
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and the Gribov parameter 7 is fixed by the gap equation: 

3Ng 2 f d A p 1 



(2tt) 4 p 4 + 7 4 L ^ 16) 

It is clear that the Gribov approach is only the first step in order to consistently treat the problem of 
zero modes and the Gribov copies in a gauge fixed Yang-Mills theory. The second step is the GZ theory 
[5j [26] , which consists in a renormalizable and local way to implement the restriction to the first Gribov 
region. In fact, Zwanziger observed that the restriction could be implemented by adding the following 
term in the action 

S G z = S YM + J 4 H(A) , (17) 
where, H (A) is the so-called horizon function, 

H(A)=g 2 J d^xdS r bc A%x)[M~ l ] ad (x,y)f dec Al. (18) 

In the Zwanziger approach, the parameter 7 is fixed by the equation 

(H(A)) = AV{N 2 - 1) , (19) 

where V is the Euclidean space volume. Notice that the Gribov form factor ()13j) coincides with the first 
order of the horizon function^: 

w^-D =a ^ A) + oiA3) - m 

It is clear that the horizon function is nonlocal, but it can be localized with the help of a suitable set of 
auxiliary fields. In order to ensure that those extra fields do not introduce extra degrees of freedom they 
are introduced in the form of a BRST quarte10: 



su>„ 



stf b = 0. 



SUJ„ 



.ab 



0. 



»■ (21) 

where (<p,tp) are a pair of complex commutating fields, while (Q,u)) are anti-commutating ones. Now, the 
local version of the GZ action is then given by: 



Q 



local 
GZ 



Sym + / d x 



,^ ac \A a b,J>c -ac\jab,,bc , 'l „ fat>c/,„at> ,^a*\ a> 



-ac \ jab, be 



abci ,^ab 



-ab\ ac 



(22) 



It is quite easy to notice that this term explicitly breaks the BRST symmetry. Then, following the 
Zwanziger steps, in order to establish a local, renormalizable and BRST-invariant theory, we define a 
most general invariant action, which possesses action (|22p as a particular physical case. Such desired 
action is then given by: 



olocal-inv 



Sym + s 



d 4 x 



,—,ac a jab , be , /~\ab r>ac,„cb , jab j-.ac-cb , r^ab j. 



s Y m + J d A x[^M ah ^ -^^+r k (^:)(^ c > ce 



i jab j-^ac,cb 



\ab r-\a< 
I ftv ^ ft 



Q 



ab 



ftV 

\f\ab T^ac -cb , / Jab jab r) ab C\ a ^ 

^^liv LJ ft "T" \ J ftv' 1 fiv ~ Wftv^ 



j-^ac, ,cb | „ facd r j-^bd d\, cl 

D ft u v +gj (D^ c )ip v 

}• 



+ J, 



ab 

ftV 



j^ac-cb | „ facai j~\bd, „d\ -cl 
D ft <Pv + 91 \P,j, c )U V 



acdi j-^bd „d\-cb 



; ftu^o ftu 



(23) 



Actually, the form factor can be calculated at all orders and the result is that such coincidence occurs in fact at all 
orders. In 17 , this equivalence is proved at third order, and in [18; it is proved at all orders. 

7 Actually, the BRST quartet is composed by two BRST doublets, which has the basic structure 

sU = V, sV = 0, 

for a generic pair of fields (U, V), guarantying the nilpotence of the BRST operator, s, i.e. s 2 = 0. 
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where, the set of external sources 
forms a BRST quartet structure, i.e. 



r jab jab f~\ab f\ab \ 
sQ an — J, 



yib 

s T ab 



ab 
fiv 

ab 



ab 



(24) 



(25) 



being (J, J) a pair of commutating sources and (Q, Q) a pair of anti-commutating ones. The last term 
is a vacuum term permitted by power-counting and it is necessary to obtain the gap equation (|19|) by 
demanding that the vacuum energy, £, is independent of 7 2 , i.e. 



d£ 

dj 2 



0. 



(26) 



The original action (|22p can be recovered from (|23p when these external sources attain their physical 
values. Namely, 



J, 



phys 



-J. 



ab 



„,2 srab x 



phys 



Q 



phys 



Q 



phys 



0. 



After perform a linear shift on the ui^ variable 



{M 



-l\ab 



d u [gf bde rt c D e Jcn-^gr c D« e c 



pbdc 7-i de „e 



one can show that 



ilocal-inv 
'GZ 



phys 



s, 



local 
GZ • 



(27) 



(28) 



(29) 



It is important to emphasize here that the Zwanziger approach described above has received some im- 
provements in recent years. In [191 120] . the model was formulated in such way that the BRST symmetry 
breaking appears as a linear breaking, while in [21] . the breaking appears as a spontaneous breaking, 
instead of an explicit on^§. 

2.2 The i-particles 

As already seen, eq. ([2]) suggests that a theory presenting Gribov-like propagators, can be rewritten in 
terms of z-particles, i.e. with propagations of "particles" with complex squared masses. Now, let us take 
a look on this concept following the lines outlined in [9]. Then, we will start with a scalar field toy model 
exhibiting a confining Gribov-type propagator: 



S 



(30) 



where 6 is a mass parameter playing an analogous role of the Gribov parameter 7. The resulting propagator 
is the desired Gribov-tipe: 



(*(p)$(-p)) 



P 



p 4 + 02 



and it can be cast in a local form exactly like in the case of the GZ model: 

S = J d 4 x X - $(-d 2 )$ + <p{-d 2 )v + V29 2 $(v? -<p)- oj(-8 2 ) 



(31) 



(32) 



The renormalization of this formulation was already proven in [22] 
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As (y?, <p) form a pair of complex field we can decouple the real part from the theory. In fact, defining 



<p = -=(U + iV) 



one can write 



S 



d 4 x 



~ <!>(-<9 2 )<I» + ^ V(-d 2 )V + 2%Q 2 <S>V - uj(-d 2 )uj + ^ U(-d 2 )U 



(33) 



(34) 



From now on, we will neglect the decoupled fields (U, cu,oj) and we will diagonalize the action above by 
introducing the new field variables: 



* = ^(A + ,) 



Thus, we have 



S 



X - \{-d 2 + 2i9 2 )\ + X - 7](-d 2 - 2i6 2 )i 1 



and the propagators in terms of the new field variables stand by 

1 



(A(p)A(-p)) 



p 2 + 2i0 2 ' 



{v(p)v(-p)) 



1 



p 2 - 2i6 2 



(35) 



(36) 



(37) 



From this expression one immediately sees that the fields A and r\ correspond to the propagation of 
unphysical modes with complex squared masses ±2i6. These are the so-called z-particles of the model. 
Notice also that, despite the imaginary terms, the action (|37p is Hermitian if we require that = r\. 

We already argued that the excitations in terms of z-particles are unphysical. Nevertheless, following 
[9], physical states can be introduced by constructing suitable composite operators out of the fields (A, 77) 
which exhibit desirable analyticity properties, as encoded in the Kallen-Lehmann spectral representation. 
Such composite operators are obtained by requiring that the i-particles fields enter in pairs, i.e. the de- 
sired operator contains as many fields of the type A as of the type rj. It ensures that in the corresponding 
correlation function only complex conjugate pairs of z-particles propagate in the Feynman diagrams. In 
the present case, the simplest example of a local composite operator with the required physical properties 
is O(x) = X(x)rj(x). In [9] it was shown that the correlation function {0(p)0{—p)) has a well defined 
Kallen-Lehmann spectral representation^: 



(0( P )0(-p)) 



d 4 q 1 1 

(2tt) 4 (p - q) 2 - 2id 2 q 2 + 2i0 2 



T + p 



with 



P(r) 



(4vr) 5 



16(9 4 



r = 26 2 



(38) 



(39) 



Furthermore, a model with interacting i-particles is dealt in [23J 



'Here we present only the 4d result. The results for the spectral density in 2d and 3d can also be found in [9]. 
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2.3 The replica model 

Let us now describe, in few words, the so-called replica model. It was first introduced in [8] as an 
alternative way to solve the Gribov problem. Here, we present the replica model by starting with a scalar 
field with a quartic self-interaction term: 

l$(-d 2 + m 2 )$ + 9 $ 4 



S = I d 4 x 



Then, we define a replica of the action above, 



S' = / d 4 x 



^&(-d 2 + m 2 )$' + 5 $' 4 



(40) 



(41) 



with the same parameters m and g. The two theories interact by a soft term depending on a free 
parameter, say fi 2 , and the replica model is then written as 



Sreplica — I d X 



- <S>(-d 2 + m 2 )$ + ^ &(-d 2 + m 2 )<S>' + i/j? + g $ 4 + g $ /4 



(42) 



The replica model (|42h enjoys a symmetry which guarantees the existence of an unique mass parameter 
m and an unique quartic coupling constant g for both sectors of the theory (the original starting point 
model and its replica). This symmetry is often called mirror symmetry and its given by: 

$ $' , $'->■$. (43) 

Notice also that, when m 2 = 0, the propagator of $ field is of the Gribov type, with [i 2 playing the role of 
the Gribov parameter! 1 "^!. Thus, it can also be diagonalized in terms of i-particles. In the deep ultraviolet 
regime, the two theories completely decouple and we the obtain that the <J? field is said deconfined 
exhibiting a Yukawa-like propagator. The fi 2 parameter might be fixed by a gap equation like in the 
GZ model (US]). 

3 Taking a closer look at the SU(3) group 
3.1 General considerations 

The SU(N) group is the group of the N x N unitary matrices with determinant equals to one: 

SU(N) := { U | UU ] = 1 , det(ET) = 1 } . (44) 

The matrices U S SU(N) can be written as 

U(u) := exp (iuj a T a ) (45) 

where ui a is a parameter, the label a runs from 1 to (iV 2 — 1), and T a are the generators of the group, 
obeying the following relations: 

[T a ,T b ] = if abc T c , 
{T a ,T b } = ^5 ab + d abc T c . (46) 

In equation (|46p [ , ] stands for the commutator, while { , } for the anti-commutator; f abc are structure 
constants, which are anti-symmetric by odd successive permutations, i.e. 

jabc jbac jacb jcba . (47) 

and d abc are the components of the completely symmetric invariant rank-3 tensor of the group. 



""For m 2 ^ the propagator behaves like the most general RGZ model [?]• 
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3.2 The SU{2) group 

Before discuss the SU(3) case it is useful to spend a few words on the 5X7(2) case. In this case there are 
three generators and they are related with the Pauli matrices, a a , as follows: 



T a = T , (a =1,2,3), 



where, 

Also, we have 
and than 



a 



1 

1 



a 



-i 

1 



a 



1 
-1 



fabc ^abc 



a" 1 a" 

y ' y 



IS 



abc 



a 



d abc = . 



a a a b 



-8 ab 
2 



(48) 

(49) 
(50) 
(51) 



3.3 The SU(3) group 

In the SU (3) group there are eight generators associated with the Gell-Mann matrices^], A a , in an analog 
way as the Pauli matrices for SU(2), 

T a = y, (a = l,...,8). 
The Lie algebra of the generators of 5C7(3) is than given by 



(52) 



A 



A a \ b 



- 5 ab + d abc — 
3 2 



(53) 



Here, the generators obey the anti-symmetric and symmetric Jacobi identities: 



[X a ,[X b ,X c ]} + [X c ,[X a ,X b }} + [X b ,[X c ,X a ]} 
[A a ,{A 6 ,A c }] + [A c ,{A a ,A b }] + [A 6 ,{A c ,A a }] 



0. 
0. 



These identities give rise to the following useful relations: 

jabc jcde _|_ jade jceb _|_ yaec jcbd 
jode^jdbc _|_ jedejdab _|_ jbdejdca 



. (54) 

In particular, the second equation of ([541 is satisfied for 5C7 (N > 3) and the first one is valid for any value 
of N. Another interesting feature of the 5L7(3) group is that its generators can be grouped in sets which 
obeys the same algebraic properties of the Pauli matrices determining then three SU(2) subalgebras. 
Defining 

A± = ^(v / 3A 8 ±A 3 ), (55) 
the three SU(2) groups embedded in the 5C7(3) are given by: 



SU (2)i 
5C7(2)„ 
5C7(2)i„ 



Ai A2 A3 

y 2 ' 2 

A4 A5 A4 
2 ' 2 ' 2 
A6 A7 A_ 

yyy 



(56) 



The Gell-Mann matrices and the structure constants of the SU(3) can be found in Appendix B. 
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The gauge field can be expanded in a base of the SU (3) generators as 



^=E4y= 2 + A i A + + 4 A - ' ( 57 ) 

0=1 a^3,8 

where 

^ = Kt! 4±4 )- (58) 

It is interesting to comment here that the pure ^4-field sector of the SU(3) Yang-Mills action, i.e. 



Afield = 7 / d 4 xF° u F° u , (59) 



is left invariant by the transformation: 



12 



(A~^ , A^, A 5 ^) — > (— 4 , 4> 4) , 

(A~,A%Al) -> (-A+,-4,-4). (60) 

Notice that these transformations interchange the groups SU(2)n and S'{7(2)in and it reminds us the 
mirror symmetry of the replica model described in the preceding section. However, we can not identify 
these two subsectors of the theory as a replica model, at least not as originally conceived in [8], because 
of the presence of hard interaction terms. Furthermore, it is important to anticipate here a result that 
we will use later. In next section we will introduce the operator d abc A b ^A c ^ coupled to an external source 
J a . A possible condensation of this operator will give rise to a non-zero expectation value in the vacuum 
of J a |^| It implies that we have to choose a direction for J a and the only direction that maintains the 
symmetry (I60j) is the direction 3. Then, choosing ( J a ) = m 2 <5 a3 we have 

i(r) ^44 = (44 - a-a-) + v£ (44 - 44) + ^ (44 - 44) . (ei) 



4 The model: The SU(3) Yang-Mills with dimension 2 condensates 

We will begin with the pure Yang-Mills action in the Landau gauge, eq.([4|). Taking this action into ac- 
count and trying to do not change the fundamental behavior in the ultraviolet regime we limit ourselves 
to operators of ultraviolet dimension 2. These type of operators give rise to soft breaking terms. In 
particular, the mass operator ^ s we ^ understood in the context of the local composite operator 

method (LCO) |24j and is responsible for a Yukawa type propagator. We will focus into another dimen- 
sion 2 operator. One that is only possible into SU(N > 3), the operator constructed with the symmetric 
structure constant and the gauge field 4> *- e - \d abc A b A c At this point it is important to emphasize 
that although we introduce the operator consistently with the study of the renormalizability via a BRST 
doublet of sources. We will be interested in bringing the sources to their physical value, following the same 
procedure of Gribov-Zwanziger |25} 126] . It is perfectly possible the introduction of these two operators in 
the quantum action and study the renormalizability of both, it is done in Appendix A. In order to turn 
clear the importance of each operator and the relation of then to the two type of solutions coming from 

12 Of course that the symmetry can be extended by the ghost and gauge fixing sectors of action Q 

13 We are not calculating the value of this condensate. We limit ourselves in this work to study the consequences of its 
existence. 
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the Schwinger-Dyson equations, the scaling type and the decoupling one, in section 3 we will discuss the 
relation of these solutions and the operators \d abc A^A^ and ^A^A®. 



In order to introduce the operator ^d^A^A^, it is also necessary to introduce a BRST doublet of 

(62) 



sources 



s\ a = iJ a , s J a = . 
Then, the BRST invariant action including the dimension 2 operator ^d abc A^A 1 ^ is given by: 



Syu + s I d*x ( - \ a d aoc A ^ 



* tcl jabc Ab ac , \a jabct p\ b\ ac 9 nabc jcde 



Sym + d^x[ -rd^A^ + X a d (d^c^A ^ - lf abc d cde AVA^ + -3 a J a 



+ a I d 4 x{ W'd^Al + X a d ll D ab c b 



(63) 



The a terms are necessary by algebraic renormalization consistency, due to the fact that there is any 
symmetry to exclude such terms 14 l. The source J a can have a nonzero expectation value. In fact, setting 
(J a ) = m 2 5 a3 , the action (163ft gives rise to a propagator of the form: 



(A-(k)Al(-k)) 



2 1 5 1 7 

_■_ 6 ai § bi + V- 5 ai 6 bi + ±_ gaigbi 

k 2 ^ k 2 + i^^ k 2 -i^^ 



U2 



with 9^ u (k) being the transverse projector, 



V3k 4 + ^- 



(5 a8 5 b3 + 5 a3 5 b8 ) 



k^k u 
6 w p~ 



(k) , (64) 



(65) 



This calculation is done in the Gell-Mann representation [TT], see Appendix B. Although explicitly break- 
ing the group, this is a soft breaking and in the ultraviolet regime is expected that the group structure 
is recovered. It is important to stress here that in spite of have two massless poles in the propagator it 
is not necessary that we have physical particles directly associated to this propagator. The interaction 
terms will mix the i-particles with the massless ones. This property can be easily observed if we rewrite 
the propagator as: 



(A a Jk)A b u (-k)) 



1 



(k 2 8 ab - im l d am ) + 



2 iab3\ 



fc 4 + ^ 



(5 a8 5 b8 + 5 a3 5 b3 ) + % 



777 



4(fe4 + !£) 



A: 2 



i=i 



3k* + ?f 



m 



1 



3^/3 k 4 + ?f 



(5 a8 5 b3 + 5 a3 5 b8 ) 



OUk) , (66) 



where d ab3 are symmetric structure constants in the 3 direction. In particular, for further use, we also 
display: 

(Al(k)AK-k)} = (Al(k)Al(-k)) = ±0^(k), 



14 Also, in Landau gauge these extra terms can be absorbed by performing a linear shift of the fields variables b a and c a 

— > o + lad , c — > c + ictA . 
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-k)) 


= (Al(k)Al(- 


-k)) 


{A%k)Al{- 


-k)) 


= (Al(k)Al(- 


-k)) 


(Al(k)Al(- 


-k)) 


= (Al(k)AU- 


-k)) 


(Al(k)Al(- 


-k)) 


= (A*(k)AU- 


-k)) 



k 2 + V- 



(k), 



k 2 _ im l 6 ^ ' 



fc 4 + 



' /it/ 



.m 



(AO 



(67) 



At this point some considerations about the explicitly breaking of the SU (3) are necessary. First, as was 
presented in section 13, 3\ the Gell-Mann matrices grouped as in (|56[) have the same algebraic properties 
as the Pauli matrices and so determine three natural SU(2) subalgebras. So taking 



hi 



A 3 , h 2 = A 8 , ei = A 1 ±iX 



A 6 ± iX 7 , e 3 ± = A 4 ± iX" 



it is easy to observe that the e\ obeys the following algebra 



[hi,h 2 
\/3/i2 - hi,e\ 



0, 



\/37i2 + /ii, ei 



Te±> 



\/3/i2 — hi,e± 



1 



1 2 

=f 5 4, 



\/3/i 2 + /ii,e 3 t 



±-e 



V3h 2 — hi,e± 



1 3 



2/ti, [e^,e?.] = Vzfa-hi, [e 3 + ,e 3 _] = V3h 2 + 



This give rise to the well know weight diagram [TT]. Moreover, from the properties of SU(2) representations 
we know that 2p = mi;\/3g — p = m 2 ;^/3q + p = WI3, where (p,q) corresponds respectively to the 
eigenvalues of hi and h 2 ordered as points in M. 2 . Assuming that the trace of the propagator must be 
real, the natural choice of the direction of the breaking in hi give rise to two goldstone bosons associated 
to e^_. The remaining 2 sets of SU(2) corresponds to a similar structure as the replica model [8J. It 
is important to stress here that this is not a model for confinement in SU(3). This is essentially an 
alternative mechanism for Gribov that presents confinement in the remaining two SU (2) groups and has 
some defined observables associated to the remaining SU(2) groups. In section 6 we will return to these 
issues in order to define the relationship between observables and the remaining group structure. 



5 Taking into account the Gribov copies. The scaling type solution 
for the gluon and ghost propagators 

In order to offer a better understanding of the model, a dynamical framework for the parameter ( J a ) = 
m 2 S a3 should be provided, i.e. we should be able to establish a gap equation for that parameter as 
in the Gribov approximation [4] or in the GZ theory [U [26] . The most immediate way to achieve a 
meaningful gap equation for m 2 is following the steps detailed in Q— > (jl6|) . also is relevant to remember 
that d aac (J c ) = Q 15 l . The no-pole condition is implemented by stating that 

a(0,A) = l, (69) 

with the propagator (164p . which yields the gap equation determining the parameter m 4 , or equivalently 
(J a )(J a ), where dimensional regularization, D = 4 — e, has been employed. Therefore, at one-loop order, 

15 This general property can be easily seen, for example, in the Gell-Mann representation for the generators. 
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the gap equation yields 



3 f d D k 1 , 1 f d D k 



m 



A 



1 9 \j {2^r k A + r^ + 9\ Q J {27t)D A;2 ( fc4 + 1 ni)( A;4 + ^)' ( 7 °) 

which enables us to express the parameter m 4 as a function of the coupling constant g. It is clear that the 
second integral is absolutely convergent and do not change the fact that m is determined as a function of a 
regularization mass A. Also from the gap equation (|69jl it follows that, for small values of the momentum 



a(k 2 ,m 4 ) 



= 1-Ck 2 , (71) 



and we recover the ghost propagator which is enhanced in the infrared {c a (k)c b (—k)} ~ -r%- We have thus 
recovered the so-called scaling solution, i.e. a suppressed gluon propagator which vanishes at the origin, 
and enhanced ghosts, which corresponds to the solution of the GZ theory. 

5.1 Taking into account ^A^A^ and the decoupling type solution. 

Recent lattice numerical simulations J27J ESI EH EH Ell E31 E2] indicates a gluon propagator which is 
suppressed in the infrared and which attains a finite non-vanishing value at zero momentum, while the 
ghost propagator turns out to be not enhanced, i.e. < c a (k)c b (—k) >« p-. This behaviour is know as 
the decoupling solution and has also been obtained from the analysis of the Schwinger-Dyson equations 
[34} l35l 136] . This solution appears in the Gribov-Zwanziger theory when the dynamics of the lacalizing 
fields is taken into account. In our model this behaviour is associated to the mass operator ^A^A^. Before 
presenting the equations that characterize the stability of this operator it is important to emphasize that 
the tad pole presented in (|104|) despite being ultraviolet convergent has problems in infrared, but not 
necessarily at k = 0. The best chance to solve this divergency is introducing a mass term in the action. 
The action £ with a mass term is given by: 

V 2 I * 



S M = S + ^-y d*xA«A«, (72) 

where £ was defined in (I63p . The BRST variation of the mass term, turns out to be proportional to the 
equation of motion of b a i.e 

S(£ M ) = 0-^jd A x (d,c a )A« = -in 2 J d 4 xc J -^ , (73) 

modifying the Slavnov- Taylor to 

<S(£„) = I d x\ — - — - H — - + ib — - + — - + iu c — -} 

K ^' J x Sn<*5A% i 5L a 5c a 5c a J 5X a * 5b a 1 

= 0. (74) 

It is now clear that the gauge propagator changes to a propagator in close relation to the one obtained 
from the refined Gribov-Zwanziger theory [371 [71 [38] . The new propagator is given by: 

< A^k)A b v (-k) > = [(— 1 -^— r )((k 2 + fW b ~ im 2 d ab3 ) 

^ )V* s ai 5 H - -( k2 + f " 2 4 )(5 a8 5 b8 + 5 a3 5 b3 ) 

4((fc 2 + M 2 ) 2 + ^) k 2 + v 2 '^ 3\k 2 +n 2 ) 2 + ^- 

+ iH^h 1 )^ § b3 + § a3 6 b 8)}]e (75) 

V33 v (A; 2 + ^ 2 ) 2 + ^ n >n M V ' 
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with m 4 defined in (|102p . This expression has a a finite nonvanishing value at zero momentum character- 
izing the decoupling type solution. It is important to emphasize here that lattice simulation results that 
have obtained the behaviour of the propagator are making use of the trace of the propagator. Remem- 
bering that d aab = 0, the dominating term that lattice is capable to see is k +fl a . It is clear that 

the dynamical origem of this mass parameter needs more explanation and we pretend to do this into a 
future work. 



6 Local composite operator and the Kallen-Lehmann spectral repre- 
sentation 

Remove the gauge fields of the physical spectrum of the theory is not enough for a model that attempts 
to offer an alternative, at least in part, to the original Gribov question. It is also necessary to present 
a candidate for physical observable that displays the Kallen-Lehmann spectral representation and cor- 
responds to an invariant composite operator. It is pointed out in [SJ [5] that a local composite operator 
constructed with i-particles [9] has one-loop correlation function that exibits the Kallen-Lehmann spectral 
representation, i.e. 

Z{P ^ ~ I (2vr) 4 ((p - k) 2 + iy/2v 2 )(k* - iV2v 2 ) ' ^ 
exhibits a spectral representation, as: 

poo 1 1 

X(p 2 )-X(0)= / dvp{v){ —3- -}, (77) 



where the spectral density p(v) = —L^. Vv' 2 - _8iA - g positive in the range of integration. This property help 
us to find an operator that has the desired analyticity properties. Now let us return to the breaking of the 
SU(3) group in order to understand the mechanism that permits the existence of physical observables. 
First of all, one of the desired aspects for one observable is not only that has Kallen-Lehmann spectral 
representation of a particle but also does not carries color index and be gauge invariant. It is clear that 
the last requirement is not possible due to the explicitly breaking of the BRST symmetry. This is also a 
problem in the original Gribov- Zwanziger due to the same problem, i.e. the explicitly symmetry breaking. 
In our case another question emerges, the explicitly breaking of the group structure. Fortunately the last 
question will be the answer in order to obtain observables that does not carry color index. At least of the 
remaining group structure. 

Analyzing the i-particles concept presented in H] we can see that it is impossible to obtain from a 
Gribov type propagator an observable that do not corresponds to an integral of the type (JT6j) . This does 
not only suggests that the observables must be constructed taking into account i-particle type correlators 
in order to obtain a real particle pole, but also indicates that a mechanism in order to do that must mix 
two different types of particles as done in [8]. The natural candidates that emerges as possible observables 
can be associated to the remaining 2 sets of C(SU (2)), in particular it is convenient to define the quantities 
E 3 + = j^( F U + %F%), E*- = -L(F^ - iF^), Ef+ = ^{F% + %F%,) and E 2 ~ = ^(F*, - iF%,) that 
obeys the same algebra as presented in (168 j) . From these operators it is convenient to define one possible 
candidate to observable as: 

= E 2 +EH + El-El- 

4> = F ^u F llu - F % F l,v ( 78 ) 

It should be noted that in spite of the breaking of the BRST symmetry, the candidate to observable must 
be a singlet. It is important to remember that the E^ v essentially obeys the same group properties as 
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defined in (|68p . Remenbering that we are performin the breaking in diretion 3 or in the notation presented 
in (|68|) hi, the algebra between h% and 4>4 is given by: 

[^,41 = ^4 [^4] = ±~ e ±' (79) 

making the algebraic result [hi, e+e+ + ele?_] = easily to obtain and proving that this operator is invari- 
ant under the remaining group charge defined by the algebra (|79j) Taking into account the propagators 
as: 



< Al(k)Al(- 


-k) > 


= <A%k)Al{- 


-k) >= 


TpOfivik) 




< 4(k)Ai(- 


-k) > 


= <A%k)Al{- 


-k) >= 


1 




k 2 + i^ 


< A%k)Al{- 


-k) > 


= <Al(k)Al(- 


-k) >= 


1 




k 2 -i^- 


< Al(k)Al(- 


-k) > 


= <A*(k)A*(- 


-k) >= 






k A + ^- 


< Al(k)Al(- 


-k) > 


= <Al(k)Al(- 


-k) >= 


m 2 1 

- % —F^ a® 

Vzk 4 + ?4- 



4 v fJ,V\ n 't, 



(k), (80) 



the first candidate which has the desired one-loop correlation function [9] is given by: 

p(v) 



< </)(k)<t)(-k) > = 12 / dv 



oo 

->■■ l' [r > , 



,2 V 



+ k 2 ' 



y/v 2 - m 4 (u 2 + m 4 ) 

= 2 • ( 81 ) 

There is another observable with different value for the mass of the condensate. In order to understand 
the second candidate to an observable, i.e. presents the Kallen-Lehmann spectral representation of a 
particle, it is convenient to consider the quadratic part of the action, in particular the term with A^ and 

4- 

Ss,s = J d'x{\Al{-d 2 )Al + \AHr&)A* + m^AjAj}, (82) 

where we have already taken into account the Landau gauge conditions, d^A^ = and d^A^ = 0. This 
sector of the action can be diagonalized trivialy making the field redefinition. 

u, = ^(4 + 4) 

V, = -J=(-4 + 4). (83) 

Therefore 

^3,8 = J d'x{hj,{-d 2 + t^=)U, + \v,{~d 2 - i^)Vj, (84) 

which describes again the i-particle structure. In order to write a physical operator that has the desired 
structure, it is sufficient to study the operator 

X = id^Uy - d v U^V v - d v V^), (85) 
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where its two-point correlation function can be cast in the form of a Kallen-Lehmann spectral represen- 
tation of a physical particle, 



Again it is convenient to observe that this operator can be associated to the abelian subgroup of the 
original SU(3) group and it is clear that the mass pole only depends on the mass gap of the i-particle. 
Again it is important to stress that due to the fact that the insertion is a soft broken, the group symmetry 
is recovered at the limit k — > oo. Here is important to comment that a dynamical symmetry breaking 
could be the answer to a BRST invariante observable. In this case the full BRST operator carries not 
only the gauge fields but also the auxiliary fields. In this case it could be possible to define a BRST and 
group invariant colorless physical observable. 

7 Conclusions 

In this work we have studied the SU(3) Yang-Mills theory in a Landau gauge with a soft mass term 
proportional to the symmetric structure constant. This soft mass term alow us to treat the Gribov problem 
in a local renormalizable action without the need of auxiliary fields like in the Gribov- Zwanziger theory. 
The confining behavior is thus induced by the soft mass term and the relation with the scaling solution 
is presented in close analogy to the Gribov mechanism and to the replica modcl[8j. Also the introduction 
of the diagonal mass term and the relation with the decoupling solution is discussed. Moreover this can 
open the possibility of a more general relation between confining behavior and the existence of gauge 
condensates in the infrared regime. In fact, it would be interesting the possibility of obtaining the value 
of the mass gap by the local composite operator method and study the relation between the extremun of 
the effective potential and the Gribov mass gap equation. 

The result we have obtained suggests further that the restriction to the first Gribov region could be 
implemented into other gauges by a simple soft mass term involving the symmetric structure constant 
that breaks the original group into a more simple group struture with a replica. Further studies aimed 
at establishing such a connection, in particular for the maximal abelian gauge, in a more precise way is 
called for and is a topic of current investigation. 

Another important point analysed is the physical operators associated to these propagator. Due to 
the structure of the propagator this is a difficult task. The simple i-particle structure is not present 
in these model and the mechanism associated to the i-particle is much more complex and involves the 
group structure itself. This is the price to be paid in order to do not double the number of gauge fields 
or introduce auxiliary localizing fields. Nevertheless we find possible candidates of physical operators 
associated to this propagator in which the group structure is present in the ultraviolet regime. It is 
argued that a dynamical breacking mechanism could restore the full BRST invariance of the observable, 
in particular for F^ v (x)F£ LV {x) in the target remaining group structure. The construction of all this 
mechanism for a target remaining group like an SU(3) replica is another topic to be investigated into 
future works. 
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A Renormalizability 



A.l Equations Compatible with the Quantum Action Principle 

The full set of equations compatible with the quantum action principle |12j is given by: 
• The Lagrange multiplier and the antighost equation. 



sr +d ^ ~ (87) 

The integrated ghost equation 
£ a (E) = A a 

A" = J d'x{ g r k (nl, + ad„\ t )Al- g r k (LV + ie\ i f)}. (88) 

Before presenting all the equations compatible with the quantum action principle it is relevant to 
note here that the term |j a j a generates a linear breaking in the ghost equation. This will gave us 
the information that this term does not renormalizes. 



Slavnov- Taylor 



, , 4 . <5£ 5E <5E <5E ,„5E n SE _ 



5(E) = 0. (89) 

• Rigid symmetry 

- Jdxgf i\ SA% +c 6cC +b §bc +c g _ c+ il^^+L 5Lc +X §xc+3 S .J 
W a (Z) = 0. (90) 

• SL(2,R) 

m - J (91) 
A. 2 Stability of the quantum action 

The next step is to characterize the most general counterterm that can be freely added to all orders 
in perturbation theory respecting all the symmetries presented previously. Following the set up of the 
Algebraic Renormalization [T2], we perturb the classical action E by adding an integrated local polynomial 
in the fields and sources, E coun i, with dimension bounded by four, and with vanishing ghost number. The 
perturbed action (E + eE cou „i), where e is an expansion parameter, fulfills, to the first order in e, the same 
Ward identities obeyed by the classical action E, i.e equations (|57]) - (|9T|) . 



Pz^count — 0, TUT - — 0) (t^T + ^77^)^ count — 0, 



S^count ~ / 

5b a 

G a ^>count = 0, W a T, count = 0, RIl coun t = 0, (92) 



5b a y 5c a 'SO,*' 
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where /3s is given by: 



d A x{ 



SE 6 



+ 



SE 5 



5E 5 SE 5 



SA° 6A° (50° 



- + T"F — ? 1" r rr + ib a — + ij a TT- } 



5c a 



5X° 



(93) 



thus, taking into account the general results on the cohomology of Yang- Mills theories, the most general 
invariant counterterm is: 



-'count 



-i 



(94) 



The integrated term corresponds to the nontrivial part of the cohomology of /3s, while A -1 is an integrated 
polynomial in the fields and sources with ultraviolet dimension 4 and ghost number -1. The ultraviolet 
dimension and ghost number of all fields and sources. 



fields and sources 


A a 


c a 


c a 




A a 


f 


% 


L a 


UV dimension 


1 





2 


2 


2 


2 


3 


4 


Ghost number 





1 


-1 





-1 





-1 


-2 



Table 1: Quantum numbers of fields and sources. 

Applying all the constraints given in (j92|) and observing the ultraviolet dimension 4 and ghost number 
— 1, we obtain for A -1 : 



A - 



J d 4 x{ ai (n« + + a 2 ad^\ a Al + ^\ a d abc AlA c ^} 



(95) 



We see thus that E coun t contains 4 free independent parameters, namely (ao, a\, a 2 , 03). These param- 
eters can be reabsorbed by means of a multiplicative renormalization of the gauge coupling constant g, 
off the parameters (a, e) and the set of fields and sources <fi = (A", c a , c a , b a , f2^, L a , A a , j a ) according to 



E(g, a, cf)) + HE 



count 



E(g ,a ,(j)o) + 0{h 2 ) , 



(96) 



with 



90 



Z g g, ao = Z a a, eo = Z E e 
1 

72 40 a _ y a _ y m -a _ y-a 



(97) 



and 



^1 



Z A 



1-fry, Z Q = 1 + - «2 - 03), Z e = 1 + h(2a - 2a 3 ), 

! + ^( a o + 2ai), = Zc = = 1 - fry , Z 6 = 1 - ^(a + 2ai), 

i-My-y-^), = 1 - h(a - a 3 ), 



or directly in terms of the multiplicative relations between the renormalization factors Z 

Z h Z\ 



- - 2 

— 1) ZcZ q Z\Z c = 1, Z £ Z; = 1. 



(98) 



(99) 



In order to clarify the importance of such relations and the Gribov type propagator, let us follow the 
Zwanziger prescription and set the sources X a ,j a respectively to 0, < j a > with < j a >^ 0. Taking this 
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into account it is clear that the action is BRST invariant up to a soft breacking term proportional to 

< f > 

sS(^,c a ,c a ,6 a ,A a = 0,L a = 0,fi£ = 0,j a =< j a >) =< f > A a 
A a = j d 4 x{d abc d^c b A^-p abc c b d cde Af l A e fl + ad ll Dfc b }, (100) 

just like in the Zwanziger procedure to the Gribov problem and due to the fact that we dont have localizing 
fields, the limit < j a >—> clearly recover the pure Yang-Mills in the deep ultraviolet region. 

A. 3 Nonrenormalization of the soft mass term, simple one loop prove 

Following closely the arguments presented by Sorella in [8] we will prove that there is no one loop correction 
to the soft mass term. The argument is based on dimensional regularization with minimal subtraction 
and the fact that, at least at one loop, the tadpole diagram in the two point A — A is related to an integral 
of the type. 

< lD X_ f adc fb ec < A d {k)A t(k) >, D = 4 - e (101) 



(27T) 



which can be rewritten as 
d D x r 



— _ — -) - if adc f bec d de3 m 2 { 



4 



+ f adc f be \ {k r + ^ Y^ 5di5ei - V-^^^ +sa ^ 
+ i $l { -^^ ){5d85e3+5d35e8))} (io2) 

Now in order to prove that all terms do not contribute to the mass, at least in first order of perturbation 
theory, let us analyse the different terms of this integral. The first term give rise to: 

d D x k 2 f d D x 1 f d D x // 

(2ir) D k 4 + ^ ~ J (2^F " J (2vr) D /fc 2 (fc 4 +/x 4 )' ( } 

corresponding respectively to an integral that is zero by dimensional and a power counting ultraviolet 
convergent one. The second term to be analysed is of the form: 

(104) 



d D x 1 f d D x 1 f d D x fT 



(2ir) D k 4 + /i 4 J (2^) D k 4 J (2tt) d k 4 (k 4 + fi 4 ) ' 

Recursively the same argument is applied to the third integral, showing thus that no divergent terms 
proportional to < > arise. Also, from the absence of one-loop counterterm of the kind < j a > 
\d abc A b A c in it follows that 

< f > ^ c (4) (^)o =< f > \d ahc A\A% (105) 

so that we obtain 

< f > = Z <3> < f >, Z Kj> Z A = 1, (106) 

meaning that the renormalization factor of the soft parameter < j a > can be expressed in terms of the 
gluon renormalization factor Za- El 



16 A purely algebraic proof, valid to all orders, of the non-renormalization properties of the soft parameter is under inves- 
tigation. Also the possibility of obtaining the value of < j a >, by the local composite operator method [24l 139] . 
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B The Gell-Mann matrices and the structure constants of SU(S) 



The Gell-Mann matrices are given by 



Ai 



o 


-1 

1 




1 





























:) 


1 






















— i 





i 







h- > 











-1 




























; ) 





1 


0/ 



A 4 = , A 5 = 



A 7 = -i , A 8 = —= 1 . (107) 



The nonzero structure constants are: 

yl23 _ ^ ^147 _ _y!56 _ y246 _ y257 _ y345 _ _y367 _ }_ y458 _ y678 _ (108) 



2 , 2 



The nonzero components of the symmetric tensor <i a6c are: 

d 118 = d 228 = ^338 = ,^888 = J_ 
^448 = d 558 = ^668 = ^778 ^ 1 



2^3' 

;146 = ^157 = _ d 247 = ^256 = ^344 = ^355 = ,^366 = ,^377 = 1 (1Q9) 
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